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Section  I.  Introduction 


The  project  entitled  “Toxicity  Evaluation  of  Engineered  Nanomaterials:  Portable  In  Vitro 
Chamber  to  Study  Realistic  Occupational  Exposure  in  Biological  Systems”  focused 
toward  developing  a  fundamental  understanding  of  the  interaction  of  nanomaterials 
(NMs)  with  biological  molecules  and  developing  the  instrumentation  and  methodology  to 
mimic  realistic  exposure  scenarios  to  predict  the  toxicological  effects  of  engineered  NMs 
in  relation  to  occupational  human  exposure.  Dr.  Saber  Hussain  was  the  overall  program 
manager  as  well  as  a  Principal  Investigator. 


Section  II.  Purpose 

The  research  gap  to  develop  realistic  ambient  NM  exposure  scenarios  in 
nanotoxicology  was  addressed.  Although  many  studies  are  being  conducted  to 
investigate  the  toxicity  of  NMs,  traditional  techniques  do  not  accurately  mimic  material 
characteristics  during  exposure,  which  can  yield  contradictory  and  inconclusive  results 
in  laboratory  NM  testing  experiments.  This  study  describes  the  methodology  necessary 
to  develop  more  realistic  NM  exposure  conditions,  in  addition  to  providing  aspects  of  in 
situ  NM  characterization  and  toxicity  kinetics  and  mechanisms.  The  final  objective  is  to 
collect  relevant  toxicological  data  to  support  accurate  data  gathering  for  high  predictive 
power  during  human  exposures  risk  assessment. 


Section  III.  Studies 

A.  Portable  In  Vitro  Chamber  to  Study  Realistic  Occupational  Exposure  in  Biological 

Systems 

Studies  conducted  by  Ms.  Christin  Grabinski,  MS 

The  main  objective  of  this  study  was  to  design  an  experimental  set-up,  incorporating 
engineered  NM  aerosolization  and  in  vitro  exposure  at  the  air-liquid  interface,  to  mimic 
realistic  respiratory  inhalation  conditions.  Traditional  exposure  testing  begins  with  a 
series  of  in  vitro  assessments  that  commonly  require  engineered  NMs  be  dispersed  in 
biological  media  prior  to  exposure.  This  approach  can  alter  the  NMs  physical  properties 
(e.g.  agglomeration)  and  chemical  surface  (e.g.  protein  coating).  There  are  many 
examples  of  agglomeration  tendency  for  NM,  but  even  the  dynamics  can  be  changed; 
for  example  mercaptohexadecanoic  acid  functionalized  gold  nanorods,  which  are  stable 
when  dispersed  in  water,  but  become  very  unstable  when  dispersed  at  the  same 
concentration  in  media  (Figure  1).  Because  many  nanoparticles  agglomerate  in 
biological  media  without  the  addition  of  a  surface  stabilizing  molecule,  we  recognized 
that  mimicking  realistic  exposure  conditions  in  the  lab  is  difficult  and  will  produce  unique 
results  from  traditional  techniques.  Consequently,  lab  model  procedures  and  systems 
for  NM  exposures  require  a  redesign  to  more  accurately  assess  human  risk. 
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Figure  1.  Gold  nanorods  dispersed  in  water  versus  media.  Plot:  UV/Vis  data  exhibited  a 
flattened  longitudinal  peak  for  gold  nanorods  dispersed  in  media  indicating  significant 
agglomeration.  Insets:  Darkfieid  images  confirm  this  phenomenon  visually.  This  figure 
was  modified  from  Grab  inski  et  al.  (201 1)  [1]. 


In  order  to  address  this  challenge,  an  instrument  set-up  was  designed  and  engineered 
for  exposing  ceiis,  grown  on  a  biocompatible  membrane  with  the  ceii  surface  side 
toward  the  environment,  to  ambient  NMs  dispersed  in  the  gas  phase.  For  our  piiot 
experiments,  a  commercially  available  horizontal  NaviCyte  chamber  (Warner 
Instruments)  was  used  as  a  closed  environment  to  expose  gas-phased,  engineered 
NMs  to  eukaryotic  ceiis,  grown  on  a  0.4  micron  porous  membrane  insert  (Corning) 
(Figure  2). 
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Porous  membrane  insert,  0.4  |im 


Figure  2.  NaviCyte  horizontal  diffusion  chamber  for  air-liquid  exposure  with  porous 
membrane  insert.  Image  of  the  chamber  was  modified  from  Warner  Instruments  Catalogue  “A 
Guide  to  Ussing  Chamber  Systems.” pp.  9-10.  Available  at 

https://www.harvardapparatus.eom/hapdfs/HAI_DOCCAT_4/Ussing_Chamber_Systems.pdf 

Different  techniques  were  used  for  dispersing  engineered  NMs  in  the  gas  phase.  For 
the  first  phase  of  the  project,  engineered  NMs  were  dispersed  using  electrospray 
techniques  (TSI,  Shoreview,  MN).  The  electrospray  operates  by  drawing  a  conductive 
buffer  containing  dispersed  engineered  NMs  into  a  capillary.  An  electric  field  is  applied, 
drawing  the  liquid  into  a  conical  jet,  which  releases  ultrafine  charged  droplets.  For  dilute 
liquid  dispersions,  each  droplet  contains  a  small  number  of  NMs,  allowing  for  uniform 
gas  phase  dispersion.  Once  droplets  containing  NMs  are  formed,  the  liquid  evaporates 
leaving  particles  in  the  gas  phase.  The  charge  on  the  NMs  is  neutralized  by  a  sealed 
ionizing  source  (Polonium-210).  The  particles  are  then  mixed  with  a  flow  of  filtered  dry 
air  (95%)  and  CO2  (5%),  which  can  then  be  directed  to  the  exposure  chamber  (Figure 
3).  Future  studies  will  investigate  additional  approaches  for  NM  aerosolization,  including 
aerosolization  of  dry  powders  and  direct  synthesis  in  the  aerosol  phase.  For  the  latter, 
metal,  carbon,  and  semiconductor  NMs  will  be  synthesized  in  the  gas  phase  using 
direct  current  atmospheric-pressure  microplasma,  which  has  been  standardized  by  Dr. 
Mohan  Sankaran  of  Case  Western  Reserve  University  [2]. 
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Specifications: 

Particle  generation  rate:  >10^  particles/cm^ 
Liquid  conductivity:  0.2  S/m  nominal 
Liquid  flow  rate:  50  to  100  nL/min 
Particle  size  range:  <3  to  >100  nm 
Initial  droplet  diameter:  150  nm 


Figure  3.  Schematic  of  Electrospray  technology. 

Schematic  was  modified  from:  Product  Information  sheet.  “Model  3480  Electrospray  Aerosol 
Generator”,  P/N  1930098  Rev  A.  Copyright  ©  2001  TSI  Incorporated  Printed  in  U.S.A.  (091401  - 
-  PDF  updated  012902).  Available  at: 

http:/ /WWW.  tsi.  com/uploadedFiles/Product_  lnformation/Uterature/Spec_  Sheets/3480,  pdf 


One  of  the  key  features  of  this  technique  is  the  requirement  for  a  conductive  dispersion. 
The  manufacturer  recommends  a  conductivity  of  2000  pS/cm  for  drawing  NMs  from  a 
liquid  dispersion  using  cone-jet  mode.  We  found  that  it  is  very  difficult  to  reach  this  level 
of  conductivity  for  stable  dispersions  of  gold  and  silver  NMs.  However,  it  is  possible  to 
generate  NM  aerosols  using  lower  conductivity  dispersions  in  silver  bullet  mode,  which 
has  also  been  shown  to  be  a  stable  approach  for  generating  aerosols  [3].  For  both  gold 
and  silver  NMs,  the  conductivity  is  high  enough  to  achieve  silver-bullet  mode.  The 
results  for  various  NM  dispersions  are  shown  in  Table  1.  Silver  (Ag)  NM  dispersions 
exhibit  a  higher  conductivity  relative  to  gold  (Au)  NMs  in  the  same  stabilizing  agent.  This 
is  due  to  the  greater  propensity  for  Ag  to  oxidize  in  solution  and  produce  Ag  ions  [4].  Au 
NMs  stabilized  with  tannic  acid  exhibit  a  much  lower  conductivity.  Also,  it  was  shown 
that  serum  proteins  could  be  added  to  increase  the  conductivity  (Figure  4).  Serum 
proteins  are  typically  used  in  cell  culture  media,  so  they  are  preferred  over  typical 
buffers,  such  as  ammonium  acetate,  which  can  affect  cytotoxicity  at  low  exposure 
levels.  In  addition  to  an  increase  in  the  conductivity,  serum  proteins  slightly  increase  the 
hydrodynamic  diameter  due  to  adsorption  onto  the  surface  of  the  NMs.  However,  this 
increase  reaches  a  maximum  value  at  around  4%  serum  (Figure  4B). 
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Table  1.  Conductivity  of  Silver  and  Gold  NM  Dispersions 


Nanomaterial 

Stabilizing  Agent 

Concentration  (mg/mL) 

Conductivity  (nS/cm) 

Ag,  60  nm 

Citrate 

0.02 

660 

Au,  60  nm 

Citrate 

0.02 

197 

Au,  60  nm 

Tannic  acid 

0.02 

120 

0.035 
0.030 
E  0.025  - 
0.020  - 


■  Gold  Nanosphere-TA 
♦  Gold  Nanosphere-CA 


0.5  1 

Percent  Serum 


1.5 


Figure  4.  Effect  of  serum  proteins  on  NM  conductivity  and  hydrodynamic  diameter.  A. 
Conductivity  increases  with  increasing  serum;  B.  Hydrodynamic  increases  with 
increasing  serum  due  to  protein  adsorption. 


In  order  to  incorporate  the  electrospray  with  the  exposure  chamber,  the  outlet  of  the 
electrospray  was  connected  to  the  chamber,  with  the  remaining  flow  directed  into  a 
HEPA  filter  (TSI).  Since  the  cells  are  sensitive  to  flow,  the  flow  rate  into  the  chamber 
was  controlled  using  a  vacuum  pump  and  mass  flow  meter  (TSI  4140).  A  heat  plate  was 
used  to  maintain  the  internal  chamber  of  the  system  at  37 °C  using  a  TC-324B  dual 
channel  controller  (Warning  Instruments).  The  heat  plate  was  custom  designed  by 
Frazier  Nash  to  incorporate  connections  for  thermocouples  and  secure  fitting  of  the 
NaviCyte  chamber  (Figure  5).  Within  the  chamber,  cells  were  grown  on  Snapwell 
membranes  (Corning)  and  fed  with  media  on  the  basal  side,  then  exposed  to  air  on  the 
apical  side  (Figure  6).  The  outlet  of  the  chamber  was  connected  to  a  HEPA  filter  (TSI)  to 
prevent  release  of  NMs  into  the  environment. 
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Figure  5.  Drawing  of  custom-designed  heat  piate  for  controlling  the  temperature  of  the 
NaviCyte  chamber. 


Temperature  Probe 


Cell  culture  membrane 
(Snapwelt  Corning) 

Well  for  basal  media 
Heat  plate 


Figure  6.  Image  showing  the  internal  features  of  the  chamber. 
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In  order  to  control  the  deposition  of  NMs  onto  the  cells  in  the  chamber,  it  was  necessary 
to  modify  the  chamber.  The  design  aiiows  for  an  eiectric  field  to  be  applied  using  an 
external  power  source  for  electrostatic  deposition  of  NMs.  Adjustment  of  the  inlet  tube  to 
create  from  0  to  1 4  cm  distance  between  the  end  of  the  tube  and  ceil  membrane  allows 
for  the  effect  of  distance  on  the  electric  fieid  and  deposition  efficiency  to  be  optimized. 
This  height  adjustment  was  facilitated  by  a  spring  mechanism  around  the  inlet  tube.  An 
electrode  was  attached  to  the  inlet  tube  and  a  conducting  materiai  piaced  beiow  the  ceii 
membrane.  The  iniet  tube  and  springs  for  adjusting  the  height  are  316  stainless  steel 
(SS)  and  completely  enclosed  within  the  system,  whiie  the  majority  of  the  chamber  was 
constructed  using  poiycarbonate  (Figure  7). 


2-Dimensional  Drawing 


3-Dimensional  Model 


Inlet  Inlet 


Outlet 

Transw'&lfansI 
support  ogp 
cells  cells 


ElectriclteidrioabtiQjifip 
from  p(dTO)OT?MpiD^su 


Dialto  (Ebattolcontrol 
inlettubelfeHi^  heigh 

Spring  8prmgtolcontr 
inlettubelfetti^  heighl 


ElectricStenariosbliQ/ffi 
forgrododg  round 


ElectroStectrode 


Prototype 

Modified  Chamber 


Adjustable  Inlet: 


Number  of 
threads  at  'O' 
on  dial 

Distance 
from  inlet  to 
cells  (cm) 

1 

0.100 

2 

0.225 

3 

0.350 

4 

0.475 

Figure  7.  Modified  Chamber.  A  drawing  and  three-dimensionai  model  were  developed, 
and  the  prototype  was  machined  based  on  this  model.  The  key  features  include  the 
adjustable  inlet  tube  and  insertion  of  electrodes,  which  will  allow  for  application  of  an 
electric  field  and  dose  control. 


NM  deposition  was  optimized  using  the  electrospray  to  draw  NMs  into  the  gas  phase 
and  eiectrostatic  deposition  (Nanometer  Aerosol  Sampler)  to  deposit  them  onto  giass 
coversiips.  Electrospray  generation  and  NM  deposition  was  investigated  for  Au  60  nm 
particies.  Greater  than  6x  more  uptake  was  observed  for  NM  deposition  using 
electrostatic  deposition,  verifying  the  appiication  of  an  eiectric  field  increases  deposition 
efficiency  (Figure  8). 
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Figure  8.  Darkfield  images  showing  Au  60  nm  particie  deposition  on  glass  coverslips 
after  30  min  generation  by  electrospray  (1.7  kV;  -066  nA).  A.  NAS  (-10,000  kV);  B. 
Chamber  (0  kV);  C.  Spectral  data  for  illuminated  spots. 


Another  set  of  experiments  were  conducted  using  Ag  NMs.  Ag  NMs  (Ted  Pella,  Inc.) 
were  characterized  before  aerosolization  and  deposition  using  transmission  electron 
microscopy  (TEM)  and  Imaged  for  primary  particle  size.  They  were  found  to  exhibit  an 
average  primary  particle  size  of  41.56  nm  (Figure  9;  Table  2).  Dynamic  light  scattering 
(DLS)  was  used  to  determine  the  hydrodynamic  diameter,  which  was  found  to  be  42.44 
nm  (Table  2). 


# 


,  •  500  nm 

Figure  9.  TEM  image  of  Ag  NMs. 


The  zeta  potential  of  the  Ag  NM  solution  was  -43.6  mV,  indicating  stable  dispersion 
(Table  2).  The  conductivity  of  the  Ag  NM  solution  was  300  |iS/cm  (Table  2). 


Table  2.  Characterization  Data  for  Ag  NMs 


Primary  Diameter 
(nm) 

Hydrodynamic 
Diameter  (nm) 

Zeta  Potentiai  (mV) 

Conductivity 

(|iS/cm) 

41 .56  ±6.76 

42.44  ±0.18 

-43.6  ±  5.2 

300 
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The  conductivity  was  sufficient  to  allow  for  steady  operation  of  the  electrospray  at  1 .45 
kV  (-40  nA)  and  1 .05  LPM  flow  (~95%  air;  5%  CO2).  Ag  NMs  (~40  nm)  were  aerosolized 
using  electrospray  in  either  an  oxygen  or  nitrogen  environment  and  deposited  onto  TEM 
grids  in  the  presence  of  an  electric  field  (-10  kV).  The  results  showed  that  aerosolization 
in  an  oxygen  environment  significantly  altered  the  morphology  (Figure  10).  The  use  of 
nitrogen  gas  for  diluting  aerosolized  NMs  does  not  alter  the  morphology. 


^  100  nm  •  100  nm 


Figure  10.  Ag  NMs  (40  nm)  stabilized  with  citrate  were  sprayed  at  a  concentration  of  20 
pg/mL  with  electrospray  settings  of  1.45  kV  (38  nA)  and  diluted  in  either  air  (A)  or 
nitrogen  (B). 


For  cytotoxicity  studies,  it  is  desired  to  achieve  even  deposition  of  NMs.  Therefore,  Ag 
NMs  (nanoComposix;  1  mg/mL)  were  aerosolized  using  electrospray  (and  diluted  in 
nitrogen),  and  the  deposition  was  investigated  using  the  Nanometer  Aerosol  Sampler 
with  maximum  DC  voltage  (~  -10  kV;  1.8-1 .9  kV/cm).  The  deposition  was  investigated 
throughout  the  TEM  grid  and  found  to  represent  equivalent  dose  without  statistically 
significance  differences  for  particles  counted  in  the  middle,  center  and  edge  capture 
regions.  A  total  of  1 0  TEM  images  were  analyzed  in  order  to  collect  the  data  displayed 
in  Figure  1 1 .  The  particles/cm^  concentrations  were  calculated  using  Imaged  software. 
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Center 


Middle 


Edge 


Location  on  grid 

Deposition  (particies/cm^) 

Center 

3.1x10®  ±6.1  xIO^ 

Middle 

2.6x10®  ±8.3  xIO^ 

Edge 

2.5x10®  ±5.9  xIO^ 

Figure  11.  Uniform  deposition  of  Ag  NMs  electrosprayed  from  a  dispersion  at  1  mg/mL 
and  dituted  in  nitrogen  +  5%  CO2.  A-C.  TEM  images  from  the  center  (A),  middle  (B),  and 
edge  (C)  of  the  grid.  D.  Table  showing  data  collected  from  10  images  representing  the 
different  parts  of  the  TEM  grid  (scale  bar  represents  500  nm);  E.  Schematic  to  explain 
words  used  for  describing  deposition  throughout  the  grid. 


Before  initiating  experiments  with  cells,  A549  immortilized  lung  epithelial  cells  were 
grown  on  a  12  mm  membrane  Insert  for  7  days  and  Imaged  to  verify  confluence.  A 
cross-section  of  the  cells  imaged  using  TEM  are  shown  in  Figure  12. 


Figure  12.  A549  cells  embedded  in  resin,  sectioned,  and  imaged  using  TEM.  The  scale 
bar  corresponds  to  2  pm. 
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NM  deposition  onto  cells  was  Investigated  using  Ag  NMs  with  diameter  of  40  and  60  nm 
from  Ted  Pella,  Inc.  (Redding,  CA).  Particles  were  delivered  to  the  modified  chamber 
containing  A549  cells.  After  exposure,  the  cells  were  washed,  and  then  trypsinized,  and 
the  amount  of  Ag  in  the  wash  solution  and  cells  was  quantified  using  Inductively 
Coupled  Plasma  Mass  Spectrometry  (ICP-MS).  The  number  of  Ag  particles  was 
calculated  based  on  the  density  of  Ag  and  estimated  volume  for  the  respective  Ag 
nanoparticle  sizes.  The  deposition  of  40  nm  Ag  particles  was  significantly  higher  than  60 
nm  Ag  particles  (Figure  13).  This  observation  was  expected  based  on  the  higher 
diffusion  rate  for  smaller  particles. 

3.0E+06  - - 1 

2.5E+06 

r.  2.0E+06 

E 

u 

^  1.5E+06 

TO  l.OE+06 
5.0E+05  1 

O.OE+00  -  , 

Ag40nm  Ag60nm 

Figure  13.  Ag  NMs  (60  nm)  stabilized  with  citrate  were  sprayed  at  a  concentration  of  20 
pg/mL  with  electrospray  settings  of  1.45  kV  (38  nA). 


The  ICP-MS  results  were  compared  to  results  obtained  using  the  Multiple  Path  Particle 
Dosimetry  (MPPD)  model  for  a  6  hour  exposure.  MPPD  calculates  deposition  & 
clearance  of  aerosols  in  all  airways  of  the  lung.  It  is  based  on  published  theoretical 
efficiencies  for  deposition  and  considers  filtration  by  the  nose  and  mouth  [5].  The  MPPD 
input  includes  spherical  Ag  particles  (40  or  60  nm),  human  airway  geometry,  and  initial 
concentration  of  0.1  mg/m^,  which  is  the  exposure  limit  for  silver  dust  set  by  the 
American  Conference  of  Governmental  Industrial  Hygienists  [6].  The  data  is  for  1  breath 
(assume  12  breaths/min).  The  results  obtained  using  the  Multiple  Path  Particle 
Deposition  model  is  shown  in  Figure  14. 
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0.120 


Figure  14.  MPPD  results  for  total  mass  per  area  as  a  function  of  airway  generation 
number  for  Ag  60  nm  and  40  nm. 


Assuming  12  breaths/minute,  the  concentration  can  be  adjusted  to  a  6  h  exposure 
(Equation  1).  These  data  can  be  adjusted  to  get  the  total  mass  deposited  (Mdep)  for  any 
exposure  time  (t),  as  shown  in  equation  1.  Also,  it  can  be  converted  to  particles  per 
surface  area  (Ndep)  using  the  particle  volume  (V)  and  density  (p),  as  shown  in  Equation 
2.  To  determine  what  the  equivalent  exposure  would  be  In  a  traditional  culture  system, 
equation  3  can  be  used.  The  deposition  efficiency  (Sdep)  is  a  function  of  NM  transport  in 
culture  media.  It  can  be  derived  using  published  theoretical  predictions  for  gravitational 
settling  and  diffusion  of  NM  agglomerates  [7]. 


PPD[  ^  ]*  [  ]*  [  ] 

/t  -  - 

p  [ — ]  = - * - 1^ — ] -  Equation  1 


rP  P  \ - 1 

p  I - ]  =  _  _  '  Equation  2 

p[  ]*  ["  ]*v[ - 

V 

P  p[_] 

—  -  Equation  3 

xp[  ]  ~  ^  P 
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Table  3.  Predicted  versus  Experimental  Results  for  Deposition. 


Ndep  (particles/cm^) 

Qdtnnlo 

1 

Predicted 

Experimental 

Ag  40  nm 

2.4x1  O^'-l. 1x1  O' 

2.6x10" 

Ag  60  nm 

6.4x10^-3.0x10'' 

5.9x10" 

The  experimental  deposition  was  lower  than  the  estimated  realistic  maximum  dose 
(Table  3).  It  is  proposed  that  higher  doses  can  be  obtained  by  generating  an  electric 
field  in  the  chamber.  The  final  outcome  of  incorporating  the  MPPD  model  and  dosimetry 
approximations  for  dosing  and  exposure  of  NMs  is  expected  to  produce  realistic  results 
for  NM  toxicity  in  vitro,  which  to  our  knowledge  have  not  yet  been  produced  in  one 
complete  study.  Results  from  conventional  dosing  at  the  air-liquid  interface  {via 
dispersion  of  NMs  in  the  culture  media)  will  be  completed  for  comparison  to  aerosol 
exposure. 

Based  on  an  extensive  review  of  literature  and  discussion  with  collaborators,  we 
discovered  that  a  critical  gap  in  the  literature  that  includes  limited  research  regarding  the 
optimal  electric  field  strength  and  frequency  for  depositing  NMs.  Therefore,  we 
purchased  a  function  generator  and  power  amplifier  where  these  instruments  can  be 
applied  to  work  together  to  adjust  the  voltage  and  frequency  coming  from  the  electrical 
outlet,  and  then  amplify  the  signal  to  the  range  necessary  for  electrostatic  deposition. 
When  voltage  is  applied  using  the  power  supply,  the  electric  field  is  achieved  between 
the  end  of  the  inlet  tube  and  the  electrode  underneath  the  cell  membrane.  The  inlet  tube 
was  fitted  with  a  fine  metal  mesh  to  allow  for  an  even  field  to  be  applied  across  the  cell 
membrane.  Due  to  the  arrangement  of  the  set-up,  there  is  negligibly  low  electrical 
current  on  the  cell  cultures.  Deposited  NMs  would  yield  a  current  measured  in  femto- 
amps  which  was  10^^  times  lower  than  currents  used  to  induce  cell  damage  through 
electroporation  [8]. 

In  order  to  investigate  the  effect  of  the  electric  field  on  cell  survival,  cells  were  grown  for 
30  minutes  in  the  chamber  with  the  electric  field  set  to  2  kV/cm.  The  Alamar  Blue  assay 
was  used  to  assess  viability  of  cells  in  the  chamber  versus  cells  grown  at  the  air  liquid 
interface  in  the  incubator.  There  was  no  negative  effect  observed  (Figure  16). 
Therefore,  the  power  supply  has  been  successfully  integrated  with  the  chamber,  and 
cell  viability  was  not  affected  by  exposure  to  the  electric  field. 
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Control  2  kV/cm 

Figure  15.  Cell  viability  data  for  cells  grown  in  the  chamber  for  30  minutes  with  a  electric 
field  of  2kV/cm  versus  cells  grown  at  the  air  liquid  interface  in  an  incubator.  Samples 
were  tested  in  triplicate.  A  student’s  t-test  revealed  that  there  is  no  statistical  difference 
between  the  control  and  test  sample. 


This  project  was  presented  during  a  poster  session  at  the  Society  of  Toxicology  meeting 
in  San  Francisco  on  13  March  2012.  The  poster  is  shown  on  the  next  page  (Figure  16). 
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Design  of  an  In  Vitro  Chamber  for  Occupational  Exposure  Assessment  of  NMs 

Christin  Grabinskj^^,  Mohan  Sankaran^  and  Saber  M.  Hussain^ 

^711  HPW/RHDJ,  Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH 
^Department  of  Chemical  Engineering,  Case  Western  Reserve  University,  Cleveland,  OH 
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Based  on  the  success  with  the  prototype  chamber,  a  scaled-up  multi-chamber  was 
designed.  A  schematic  of  the  scaled  up  chamber  Is  shown  in  figure  1 7. 


Model  Images  provided  by  Fi  azer-Nash  USA  modeled  in  Inventor 


Figure  1 7.  Multi-chamber  modei 
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B.  Evaluation  of  Basic  Toxicity  of  Various  Metal  Nanoparticles  for  QSAR  Modeling 

Studies  conducted  by  Nicoie  Schaeubiin,  MS  in  coilaboration  with  Jerzy 
Leszczynski,  PhD 


Quantitative  Structure  Activity  Relationship  (QSAR)  methods  are  based  on  the 
assumption  that  the  variance  in  a  given  physico-chemicai  or  biologicai  (e.g.,  toxicity) 
property  in  a  set  of  compounds  (so-called  endpoint)  was  determined  by  the  variance  in 
their  molecuiar  structures,  encoded  by  so-caiied  descriptors.  This  was  particulariy 
important  in  reiation  to  the  ethicai  issues,  costs  and  the  expected  increase  in  number  of 
engineered  NMs  that  will  be  utilized  in  the  future.  The  QSAR  method  can  be  used  to 
predict  the  toxicity  of  various  metai  nanoparticles  by  using  molecular  descriptors  along 
with  resuits  from  standard  toxicity  evaluations.  However,  in  order  to  validate  the  model, 
the  mathematical  equation  requires  data  from  many  different  nanoparticies. 

The  first  step  of  this  study  invoives  predictions  of  physicochemicai  properties  of  the 
studied  NMs,  in  order  to  deveiop  a  set  of  the  descriptors  to  be  used  to  study  their 
toxicity.  For  the  initial  investigations  five  compounds-  nickel  oxide,  silicon  dioxide,  tin 
oxide  and  titanium  dioxide-  were  selected  from  the  group  of  metal  oxides  tested 
experimentaliy.  Based  on  the  experimentai  crystai  iattice  parameters  we  constructed 
and  studied  metal  (semimetal)  oxides  (MeOx)  moiecular  ciusters.  For  all  these 
compounds  geometric  optimizations  were  done. 

Various  size  ciusters  were  considered,  i.e.,  metal  oxide  ciusters  of  size  range  from  2  to 
500  atoms.  In  order  to  obtain  cubic  clusters  having  an  appropriate  size,  the  lattice 
parameters  were  subsequently  increased  in  respect  to  all  three  dimensions.  At  first,  we 
increased  twice  one  of  the  dimensions  of  the  initiai  unit  cell.  Thereafter,  we  used  the 
initial  cluster  structure  obtained  in  the  previous  step  and  additionaiiy  increased  twice  the 
second  dimension.  Then,  we  doubled  the  third  of  the  unit  cell  dimensions,  increasing 
the  cluster  size  more  than  twice  the  initial  size.  We  repeated  this  procedure  until  we 
obtain  clusters  with  the  proper  size.  The  results  are  presented  in  Tabie  4,  where  the 
saturation  point  refers  the  size  of  the  considered  metal  oxide  nano-cluster  (in  A)  at 
which  the  eiectronic  property  vaiue  (Band  GAP  energy)  of  that  cluster  converges  to  the 
value  of  bulk  size  metai  oxide.  The  Band  GAP  energy  generaiiy  refers  to  the  energy 
difference  (in  eiectron  volts)  between  the  top  of  the  valence  band  and  the  bottom  of  the 
conduction  band  (for  example,  in  insulators  and  semiconductors). 

Having  both,  high-quality  experimentai  data  and  moiecuiar  descriptors:  (i)  image 
descriptors  (based  on  images  taken  from  transmission  electron  microscopy)  and  (ii) 
quantum-mechanicai  descriptors  (based  on  quantum-chemical  calculations),  a 
computational  study  was  carried  out.  The  first  step  of  this  study  was  aimed  at 
investigating  how  the  method  of  splitting  data  into  a  training  set  and  a  test  set  influences 
the  external  predictivity  of  quantitative  structure-activity  relationships  (nano-QSAR) 
modeis. 
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Table  4.  Experimental  Compounds  Selected  From  the  Group  of  Metal  Oxides  Based  On 
Their  Experimental  Crystal  Lattice  Parameters 


Metal  /  Semi- 
Metal  Oxide 

Saturation  point 
[A] 

Number  of 
atoms 

Totai  energy 
[eV] 

Band  GAP  energy 
[eV] 

AisOs 

12 

160 

-31466 

-4.59 

ZnO 

11 

144 

-23157 

-3.87 

Sn02 

13 

144 

-41962 

-3.85 

Si02 

10 

96 

-21060 

-3.23 

Ti02 

13 

144 

-31518 

-7.47 

For  the  initial  computational  investigations  five  alternative  splitting  algorithms  were 
selected,  namely  (i)  a  commonly  used  algorithm  (‘X:1’),  in  which  the  compounds  were 
sorted  according  to  the  ascending  response  values,  y,  and  every  xth  (e.g.  third)  was 
selected  into  the  test  set;  (ii-iv)  three  variations  of  the  Kennard-Stone  algorithm  [9]  and 
(v)  the  duplex  algorithm  [10].  In  the  next  step,  the  experimental  data  and  molecular 
descriptors  were  used  to  calibrate  and  validate  the  nano-QSAR  models  according  to  the 
golden  standards  established  by  the  Organization  for  Economic  Co-operation  and 
Development  (OECD)  [11].  We  employed  a  wide  range  of  different  regression 
techniques,  such  as  Multiple  Linear  Regression  (MLR),  Ordinary  Least  Squares  (OLS), 
Partial  Least  Squares  (PLS),  and  Principal  Component  Regression  (PCR).  In  order  to 
select  the  optimal  combination  of  the  molecular  descriptors  to  be  utilized  in  the  final 
nano-QSAR  model,  we  employed  the  Holland’s  genetic  algorithm.  The  algorithm 
minimizes  the  prediction  error  by  searching  for  the  most  optimal  combination  of  the 
descriptors. 

Individual  nano-QSAR  models  were  characterized  by  low  values  of  the  statistics: 
goodness-of-fit  (0.51  <  R^  <  0.65)  and  external  predictivity  (0.34  <  Ext  <  0.53).  The 
absolute  values  of  root  mean  square  error  of  prediction  (RMSEP)  were  ranged  from 
0.45  to  0.58  logarithmic  units  of  LD50.  The  descriptors-to-compounds  ratio  (5:1)  was 
satisfactory.  Various  statistical  parameters  for  ten  models  are  presented  in  the  Table  5. 


Table  5.  Statistical  results  for  Individual  nano-QSAR  models 


Model 

1 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

o2 

0.67 

0.52 

0.49 

0.61 

0.59 

0.52 

0.53 

0.61 

0.59 

0.52 

X 

LU 

CM 

C 

0.53 

0.41 

0.42 

0.49 

0.51 

0.41 

0.48 

0.48 

0.34 

0.41 

R2 

0.51 

0.57 

0.51 

0.53 

0.56 

0.55 

0.52 

0.55 

0.51 

0.65 

RMSEc 

0.21 

0.27 

0.38 

0.34 

0.33 

0.39 

0.31 

0.34 

0.43 

0.29 

RMSEcv 

0.44 

0.48 

0.52 

0.41 

0.59 

0.47 

0.46 

0.37 

0.49 

0.47 

RMSEp 

0.49 

0.56 

0.54 

0.58 

0.45 

0.52 

0.49 

0.58 

0.55 

0.52 

Particle  size  is  not  the  only  possible  factor  influencing  toxicity  of  engineered  NMs.  The 
following  features  should  be  also  considered:  size  distribution,  agglomeration  state, 
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shape,  porosity,  surface  area,  chemical  composition,  structure-dependent  electronic 
configuration,  surface  chemistry,  surface  charge,  and  crystal  structure.  Thus,  novel 
“nano-descriptors”  should  reflect  these  characteristics.  In  order  to  find  the  best  structural 
parameters  reflecting  the  essential  properties  of  the  studied  NMs,  we  proposed  a  set  of 
image  descriptors  (based  on  images  taken  from  TEM)  and  quantum-mechanical 
descriptors  (based  on  quantum-chemical  calculations). 

Based  on  the  experimental  crystal  lattice  parameters,  we  constructed  metal  (semimetal) 
oxides  molecular  clusters  with  the  XP  molecular  modeling  software  package.  We 
utilized  those  structures  as  inputs  for  quantum-mechanical  calculations  which  included 
two  stages:  (i)  optimization  of  the  molecular  geometry  with  respect  to  the  energy 
gradient  and  (ii)  calculation  of  the  descriptors  based  on  the  optimized  geometry.  The 
calculations  were  performed  at  the  semi-empirical  level  of  the  theory  with  use  of  PM6 
method  in  the  MOPAC  2009  software  package.  We  calculated  25  molecular  descriptors 
(Table  6). 
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Table  6.  Definition  of  quantum-mechanical  descriptors 


Symbol 

HOF 

Hof_k 

U 

TE 

EE 

Core 

SAS 

HOMO 

LUMO 

GAP 

n 

S 

CO 

EV 

EC 

O- 

Q+ 

Ahof 

Ad 

EN 

Hard 

Shift 

HOFc 

TEc 


Definition  of  quantum-mechanical  descriptors 

Standard  heat  of  formation  of  metal  oxide  nanocluster 

Standard  heat  of  formation  of  cation 

Lattice  energy 

Total  energy 

Electronic  energy 

Core-core  repulsion  energy 

Solvent  accessible  surface 

Energy  of  the  Highest  Occupied  Molecular  Orbital 

Energy  of  the  Lowest  Unoccupied  Molecular  Orbital 

HOMO-LUMO  energy  gap 

Chemical  hardness 

Total  softness 

Electrophilicity  index 

Electronic  chemical  potential 

Valance  band 

Conduction  band 

Lowest  negative  Mulliken’s  partial  charge  on  the  molecule 
Highest  positive  partial  charge  on  the  molecule 
Polarizability  derived  from  the  heat  of  formation 
Polarizability  derived  from  the  dipole  moment 
Mulliken’s  electronegativity 
Parr  and  Pople’s  absolute  hardness 
Schuurmann  MO  shift  alpha 

Standard  heat  of  formation  in  a  solution  represented  by  the  Conductor-like 
Screening  Model  (COSMO) 

Total  energy  in  a  solution  represented  by  the  Conductor-like  Screening 
Model  fCOMSOf _ 


Moreover,  based  on  the  TEM  microscopic  images,  we  have  converted  the  TEM  pictures 
to  a  numerical  matrix  in  which  each  numerical  value  corresponds  to  a  single  pixel  of  the 
original  picture  and  calculated  a  set  of  image  descriptors  reflecting  the  size  distribution, 
agglomeration  state,  shape,  porosity,  and  surface  area  for  all  studied  nanometer-sized 
metal  and  semimetal  oxides.  We  calculated  the  12  image  descriptors  listed  in  Table  7. 


Table  7.  Definition  of  Size  Descriptors 
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Symbol 

Definition  of  size  descriptors 

A 

Area 

V 

Volume 

dS 

Surface  diameter 

dV/m 

Volume/mass  diameter 

dSauter 

Volume/surface  diameter 

AR_x 

Aspect  ratio  X 

AR_y 

Aspect  ratio  Y 

PX 

Porosity  X 

PY 

Porosity  Y 

P 

Total  porosity 

Sphericity 

fcirc 

Circularity 

Having  high-quality  experimental  data  and  molecular  descriptors  that  included  (i)  image 
descriptors  (based  on  images  taken  from  transmission  electron  microscopy)  and  (ii) 
quantum-mechanical  descriptors  (based  on  quantum-chemical  calculations),  we 
developed  a  QSAR  model  following  the  golden  standards  and  recommendations  of  the 
Organization  for  Economic  Cooperation  and  Development  (OECD)  [11].  We  employed 
the  Multiple  Linear  Regression  (MLR)  as  the  chemometric  method  of  modeling.  In  order 
to  select  the  optimal  combination  of  the  molecular  descriptors  to  be  utilized  in  the  final 
QSAR  model,  we  employed  the  Holland’s  genetic  algorithm  (GA)  [12].  The  algorithm 
minimizes  the  prediction  error  by  searching  for  the  most  optimal  combination  of  the 
descriptors.  Genetic  Algorithm-Multiple  Linear  Regression  calculations  were  performed 
with  MATLAB  7.6  and  PLS  Toolbox  5.2. 

We  obtained  a  QSAR  equation  utilizing  only  two  descriptors  to  successfully  predict  the 
cytotoxicity  of  the  MeOx  nanoparticles. 


log(1/EC50)'^^’  =  2.442  -h  0.192*HoF  -h  0.429*EN 


Table  8.  Results  from  QSAR  Model. 
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ID 

Heat  of 
formation 
[HoF] 

Mulliken 

electronegativity 

[EN] 

Observed  log 
(1/LD50) 
[molar] 

S 

(3:1) 

Predicted 

log(1/LD50) 

[molar] 

Leverages 

Residuals 

ZnO 

-449.38 

8.329 

3.32 

1 

3.28 

0.437 

0.04 

1 11203 

-52.07 

6.784 

2.92 

1 

3.01 

0.278 

-0.09 

La203 

-157.72 

6.455 

2.87 

2.88 

0.196 

-0.01 

CoO 

-1675.45 

8.251 

2.83 

1 

2.85 

0.497 

-0.02 

Sn02 

-145.79 

4.685 

2.67 

1 

2.43 

0.095 

0.24 

M 11203 

-89.59 

4.998 

2.64 

1 

2.53 

0.107 

0.11 

WO3 

-728.13 

6.728 

2.56 

2.77 

0.143 

-0.21 

Bi203 

-148.46 

5.343 

2.50 

1 

2.60 

0.107 

-0.10 

NiO 

68.02 

4.472 

2.49 

1 

2.44 

0.136 

0.05 

Gd203 

-1784.98 

5.909 

2.41 

1 

2.20 

0.508 

0.21 

Sb203 

-206.73 

4.459 

2.31 

2.35 

0.091 

-0.04 

Cr203 

-235.25 

4.356 

2.30 

1 

2.31 

0.092 

-0.01 

V2O3 

-139.54 

3.239 

2.24 

1 

2.05 

0.175 

0.19 

Y2O3 

-135.28 

3.348 

2.21 

1 

2.08 

0.164 

0.13 

Si02 

-618.26 

3.813 

2.12 

2.04 

0.162 

0.08 

Fe203 

-416.51 

4.573 

2.05 

1 

2.31 

0.081 

-0.26 

AI2O3 

-599.96 

3.436 

1.85 

1 

1.95 

0.204 

-0.10 

TiOs 

-651 .58 

4.334 

1.76 

1 

2.17 

0.120 

-0.41 

The  descriptor  HoF  represents  the  density  of  the  MeOx  nanoparticles  surface  and  the 
avaiiabiiity  of  the  metal  cations.  In  other  words,  HoF  descriptor  was  reiated  to  how  much 
energy  should  be  delivered  to  detach  the  metai  and  transform  an  atom  into  its  fully 
ionized  form.  The  consequence  is  to  initiate  changes  in  ceii  viabiiity,  generation  of 
reactive  oxygen  species  (ROS)  and  the  response  of  human  keratinocyte  (HaCaT)  ceiis. 
The  second  descriptor,  EN,  was  ciosely  reiated  to  the  redox  potential  of  the  studied 
structures.  Those  oxides  having  high  EN  usually  exhibit  the  strong  reductive  properties 
(easiiy  detach  the  eiectron),  which  enhances  their  cytotoxicity  to  the  HaCaT  ceil-line. 

An  integrai  part  of  QSAR  modeiing  was  to  appropriately  describe  the  borders  of  the 
optimum  prediction  space  of  the  model.  We  verified  the  appiicabiiity  domain  by  using 
the  Williams  plot,  which  was  the  plot  of  the  leverage  values  versus  cross-validated 
standardized  residuais  (Figure  18). 
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Figure  18.  Williams  plot  describing  applicability  domains  of  GA-MLR  model.  Solid  lines 
represent  the  residual  threshold  (0  ±  3  standard  deviation  units),  and  dotted  line 
represents  the  critical  leverage  value  (h*). 


The  QSAR  model  was  characterized  by  satisfactory  goodness-of-fit,  robustness  and 
external  predictive  performance  (the  statistical  measures  are  summarized  in  Table  9). 


Table  9.  Statistical  Measures  for  QSAR  Model. 


Feature 

Measure 

nanoQSAR  model 

R2 

0.82 

Goodness-of-fit 

RMSEC 

0.175 

Q2CV 

0.70 

Robustness 

RMSECV 

0.222 

Q2Ext 

0.84 

Predictivity 

RMSEP 

0.112 

A  visual  correlation  between  the  experimental  and  predicted  values  of  log(1/EC50)  is 
presented  in  Figure  19. 
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Figure  19.  A  visual  correlation  between  the  experimental  and  predicted  values  of 
log(1/EC50) 


This  project  was  presented  during  a  poster  session  at  the  Society  of  Toxicology  meeting 
in  San  Francisco  on  13  March  2012.  The  poster  is  shown  on  the  next  page  (Figure  20). 
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Experimental  Results 


Cellular  MabtHty 


Theoretical  Model 


nanoQSAR  Model 


^og(i/EC.o)  =  2442  +  o.192*HoF  +  o429*EN 
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Table  2.  Medianical  Descriptors  Used  to  Calculate  nanoQSAR  model 


Summary 


^^uman  keratinoc)'te5  were  used  as  a  model  for  dermal  toxicity 

exposure  to  21  different  Me  ONPs. 

*The  LDje  molecular  descriptors  we-e  calculated  fcr  each 
pardclefromthecell  \iabiLity  assay. 

•The  LDjb  molecular  descriptors  were  used  to  sucoess&iEy 
develop  a  nanoQSAR  model  to  accurate])'  predict  the  toxicity  of 
MeONPs. 
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Figure  20.  Poster  presented  at  the  SOT  conference  on  13  March  2012 
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Additional  data  were  generated  to  make  the  model  more  reliable,  including  the  LSC 
(least  significant  concentration  of  particles  that  induce  ROS)  and  LD200  (lowest  dose  that 
generates  50%  more  ROS  than  the  control).  Therefore,  all  21  metal  oxide  nanoparticles 
(MeONPs)  were  assessed  to  see  if  they  generate  reactive  oxygen  species  (ROS). 
Each  MeONP  was  tested  at  3  different  concentrations  based  on  the  LD50  values 
calculated  from  the  viability  assay  in  order  to  streamline  the  process.  However, 
because  previous  results  failed  to  show  a  significant  increase  in  ROS  after  24h,  the 
cells  were  exposed  to  the  MeONPs  for  4h.  We  shortened  the  time-point  because  ROS 
is  a  fairly  early  cellular  response  and  the  cells’  antioxidant  defenses  may  have  already 
eliminated  any  ROS  generated  before  measured  at  24h.  The  results  from  the  4h 
exposure  demonstrated  significant  amounts  of  ROS  after  exposure  to  all  MeONPs, 
suggesting  that  the  cell  was  scavenging  the  ROS  produced  and  ameliorating  the  toxicity 
of  the  NPs  (Figure  21).  Further  experiments  must  be  conducted  to  calculate  the  LD200 
and  LSC. 
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^  Control  ^  Zinc  □  Zirconium 
150-  □  Cobalt  II  □  Cobalt  III  □  Manganese 


Control  CZI  Antimony  Indium  CZI  Lanthanum  ^1  Nickel 
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Figure  21.  Generation  of  ROS  after  4h  exposure  to  MeONPs.  There  was  a  significant 
increase  in  ROS  after  exposure  to  ait  NPs. 
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C.  Fundamental  Understanding  of  Biological  Compatibility  of  Nanomaterials 

Studies  conducted  by  Ravindra  Pandey,  PhD 

In  the  first  study,  the  interaction  of  protein  with  nanoparticle  was  investigated  via 
examining  the  interaction  of  tryptophan  (Trp)  with  gold  (Au)  nanoparticles.  The  results 
showed  that  the  interaction  strength  is  site-selective  and  is  dominated  by  tt  interactions 
together  with  small  ionic  interactions.  In  the  Trp-Au  bio-conjugated  complex,  the  nature 
of  Trp’s  molecular  orbitals  is  modified  by  the  hybrid  states  associated  with  Au  as  low 
lying  excited  states.  The  non-radiative  energy  transfer  to  these  states  appeared  to 
cause  the  experimentally  observed  fluorescence  quenching  Trp  interacting  with  gold. 

For  the  second  study,  we  performed  a  simulation  of  lipid  bilayers  using  first  principles 
methods.  The  calculated  results  showed,  for  the  first  time,  the  accuracy  and  reliability  of 
such  methods  in  calculating  the  structural  properties  for  DMTAP,  DMPC  and  DMPE 
bilayer  systems.  We  also  examined  the  solvation  effects,  and  find  that  the  presence  of 
the  water  molecules  is  essential  in  stabilizing  the  configuration  of  the  DMPE  bilayer.  It 
also  facilitated  a  significantly  large  change  in  its  dipole  moment  along  the  y  direction. 
The  results  therefore  predicted  a  relatively  high  activity  of  the  hydrated  lipid  layer  of 
DMPE  relative  to  that  of  DMPC  for  nanoparticles. 


1.  Protein-Nanoparticle  Interactions:  Tryptophan  with  Gold 
(J.  Physical  Chemistry  C  1 15,  22818  (2011)) 

Introduction 

The  miniaturization  of  structures  in  the  guise  of  nanoscale  materials  have  created  a  so- 
called  nanoplatform  for  applied  biomedical  research  in  developing  novel  and  improved 
diagnostics  and  therapeutics  [14].  Nanoparticles  possess  large  surface  areas  for 
biomolecular  reactions,  thus  providing  a  promising  pathway  for  combining  biomolecular 
functionality  with  unique  electronic  and  optical  attributes  to  create  hybrid  nano 
scaffolds.  Such  scaffolds  with  specific  features  are  of  particular  interest  in  recognition 
and  modulation  of  biomacromolecules,  protein  interactions,  protein-nucleic  acid 
interactions,  enzyme  activity,  and  biomimetic  reactions  [14-16].  Due  to  intriguing 
optoelectronic  properties  and  biocompatibility,  gold  nanoparticles  have  been  widely 
exploited  for  diagnostics  and  therapeutic  applications  [15].  The  long  term  fascination 
with  gold  was  substantially  facilitated  by  the  ease  of  their  size-controlled  synthesis  and 
functionalization  with  multifunctional  ligands  such  as  amino  acids,  proteins,  peptides 
and  DNA  [15].  Serum  proteins  associate  with  the  nanoparticles  in  biological  medium  to 
form  a  protein  ‘corona’  which  defines  the  biological  identity  of  the  particles  within  the 
biological  system,  thus  making  the  protein-nanoparticles  interaction  of  particular 
interest  [17].  The  adsorption  of  proteins  on  nanoparticles  can  alter  their  three 
dimensional  structure,  activity  and  biological  responses,  leading  to  unpredicted 
perturbations  in  the  system. 

The  protein-nanoparticle  interactions  and  their  effects  on  proteins  structure  and  activity 
are  generally  determined  by  various  spectroscopic  and  thermodynamic  methods.  One 
method,  intrinsic  fluorescence  quenching  is  a  powerful  spectroscopic  tool  for 
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investigating  protein  structure  and  binding  interactions  with  nanoparticles  [18-19].  The 
fluorescence  of  protein  is  due  to  aromatic  amino  acid  residues,  among  which, 
tryptophan  (Trp)  has  a  strong  fluorescence  and  high  quantum  yield  attributed  to  its  tt-tt* 
transitions  of  indole  functional  group  [20-21].  Trp  is  an  essential  amino  acid  in  the 
human-diet  and  is  known  to  form  hydrogen  bonds  together  with  weak  ionic  bonds  in  a 
conjugated  system.  Also,  the  indole  group  contains  a  benzene  ring  fused  with  nitrogen 
containing  a  pyrrole  ring,  resulting  in  a  delocalized  tt  electron  cloud.  At  the  excitation 
wavelength  of  280  nm,  the  fluorescence  of  protein  is  dominated  by  Trp  along  with  minor 
contribution  from  tyrosine  and  phenylalanine.  The  intensity  and  wavelength  of  Trp 
fluorescence  are  sensitive  to  the  environment,  and  hence  indicative  of  the  protein 
conformational  alterations  upon  binding  to  nanoparticles. 

Earlier  studies  suggested  that  gold  efficiently  quenches  the  fluorescence  of  Trp9  and 
the  efficiency  of  the  quenching  represents  the  binding  affinity  and  accessibility  of  gold 
particles  to  the  Trp  residue  in  protein.  In  the  experimental  study  of  bioconjugated  gold 
nanoparticles  with  BSA  protein,  quenching  of  the  static  mode  of  intrinsic  Trp 
fluorescence  revealed  the  formation  of  the  ground  state  complex  consisting  of  protein 
and  nanoparticles  [18,23-24].  The  moderate  binding  observed  for  gold-BSA  conjugation 
was  consistent  with  week  van-der  Waals  and  hydrophobic  interactions  between 
nanoparticle  and  biomolecules  [18,23-24]. 


We  next  investigated  the  Trp  fluorescence  quenching  by  gold  nanoparticles  using 
quantum  mechanical  methods.  We  employed  first-principles  methods  based  on  density 
functional  theory  (DFT)  to  determine  the  details  of  the  interaction  between  gold  (Au) 
nanoparticles  and  amino  acid  tryptophan  (Trp).  We  determined  the  (i)  preferred 
interaction  site  of  Trp,  (ii)  chemical  bonding  of  Au-Trp  bioconjugates  and  (iii) 
modifications  in  Trp’s  electronic  states  that  may  be  responsible  for  Its  fluorescence 
quenching  in  the  bioconjugate  configuration. 

Au32,  a  cage-like,  hollow,  gold  cluster,  8.9  A  In  diameter  was  used  to  represent  gold 
nanoparticles.  A  recent  report  on  magic  numbers  in  gold  clusters  by  using  the  shell  from 
the  Jellium  model  suggest  that  AU32  is  highly  stable  and  a  good  potential  candidate  for 
studying  cage-llke  structures  with  an  energy  gap  comparable  to  other  cage-like 
structures:  Ceo,  Auie,  AU42  [25].  The  binding  sites  offered  by  Trp  are  salt-bridge 
interaction  involving  the  carboxylic  group  (-COO-),  charge-solvent  Interaction  involving 
the  indole  ring,  and  amine  group  interaction  involving  -NH2  group  [26].  The  equilibrium 
configurations  with  corresponding  interaction  energies,  density  of  states,  molecular 
orbitals,  and  charge  density  distribution  of  the  Trp-Au32  complex  were  obtained  through 
DFT-based  quantum  mechanical  calculations;  all  of  which  help  define  the  interactions 
between  Au  and  Trp. 


Computational  Method 

All  calculations  were  carried  out  within  the  framework  of  density  functional  theory  (DFT) 
as  Implemented  In  SIESTA  code  using  exchange-correlation  function  proposed  by 
Perdew,  Burke  and  Ernzerhof  (PBE)  [27-28].  Only  the  valence  electrons  of  each 
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chemical  species  in  the  Trp  and  Au  nanocluster  were  treated  explicitly  and  their 
interactions  with  ionic  cores  were  described  by  ultra-soft  norm-conserving 
pseudopotentials  using  the  Trouiller  and  Martins  scheme  [29].  The  valence 
configurations  of  H  (1s^),  C  (2s^  2p^),  N  (2s^  2p®),  O  (2s^  2p'‘),  and  Au  (5d^°  6s^)  were 
represented  by  double-zeta  basis  sets  with  polarization  functions  (DZP)  [30].  The 
numerical  integrals  were  performed  on  a  real  space  grid  with  an  equivalent  energy 
cutoff  of  3400  eV.  The  calculations  were  considered  to  be  converged  when  the  force  on 
each  ion  was  less  than  0.001  eV/A.  Also,  total  energy  convergence  criterion  of  1 0  ®  eV 
was  used. 


Results  and  Discussion 

The  32  atom  cluster  of  Au  was  found  to  favor  a  cage-llke  structure  [31].  The  calculated 
ground  state  configuration  of  AU32  has  icosahedral  I/,  symmetry  with  hollow  shape  of  the 
cluster  as  shown  in  Figure  22a.  The  geometry  was  obtained  by  placing  one  atom  on 
each  pentagon  of  a  dodecahedron.  The  cluster  evidenced  only  two  different  kinds  of 
sites  on  the  cage:  12  sites  with  a  coordination  number  five  and  20  sites  with  a 
coordination  number  six;  i.e.  each  atom  binds  to  either  five  or  six  neighboring  gold 
atoms.  The  symmetry  Is  similar  to  the  truncated  icosahedron  Ceo,  with  interchanged 
vertices  and  planes.  Therefore,  two  bond  lengths  of  2.72  and  2.8  A  were  obtained  for 
the  five-coordinated  sites  and  six-coordinated  sites,  respectively  (Figure  22a).  The  total 
charge  density  of  AU32  is  also  shown  in  Figure  22b.  The  calculated  structural  properties 
of  AU32  are  in  good  agreement  with  the  previously  reported  values  using  density 
functional  theory  [31].  The  gap  between  the  highest  occupied  molecular  orbital  (HOMO) 
and  the  lowest  unoccupied  molecular  orbital  (LUMO)  is  calculated  to  be  1 .31  eV. 

The  presence  of  an  indole  functional  group  together  with  an  amine  and  a  carboxyl 
functional  groups  is  the  discriminating  structural  characteristic  of  the  aromatic  amino 
acid  residue  tryptophan  (Trp).  Figure  22  shows  the  optimized  ground  state  configuration 
and  total  charge  density  of  Trp  molecule.  The  indole  ring  is  oriented  perpendicular  to  the 
carboxyl  and  amine  groups.  In  the  optimized  configuration,  the  indole  ring  is  113.77- 
out-of-plane  with  respect  to  the  plane  containing  -COOH  and  -NH2  functional  groups 
i.e.  C1-C2-C4.  The  C-C  and  C=C  bond  lengths  in  the  hexagonal  benzene  ring  of 
indole  are  1 .42  and  1 .40  A,  respectively.  The  C-C  and  C=C  bond  lengths  of  the 
pentagonal  pyrrole  ring  of  indole  are  1.45  and  1.39  A,  respectively.  The  C-C  bond 
distinguishing  the  plane  of  the  indole  ring  and  the  -COOH/-NH2  functional  groups  is 
1.56  A  making  an  angle  (C2-C4-C9)  of  105.45-  with  -COOH  group  directed  towards 
the  indole  ring.  The  O-C-0  angle  of  -COOH  group  Is  1 23.23-.  The  calculated 
HOMO-LUMO  gap  and  dipole  moments  are  3.66  eV  and  4.01  debye  (D),  respectively. 
The  calculated  values  of  bond  lengths  and  bond  angles  of  Trp  are  in  excellent 
agreement  with  the  previously  reported  values  also  obtained  using  DFT  calculations 
[32]. 
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Figure  22.  The  calculated  ground  state  configurations  and  total  charge  density  of  (a) 
(AU32),  and  (b)  Trp.  The  AU32  interaction  sites  for  tryptophan  are  also  shown.  The  charge 
density  is  shown  at  1/3rd  of  its  maximum  isosurface  value. 


The  optimized  isomeric  configurations  of  Trp-Au32  conjugates  involving  (i)  -COOH  (ii) 
-indole,  and  (iii)  -NH2  functional  groups  as  preferred  interaction  sites  are  presented  in 
Figure  23.  No  geometrical  constraints  were  imposed  to  obtain  the  optimal  structures  of 
these  isomeric  configurations  of  Trp-Au32-  The  binding  energy,  bond-lengths,  dipole 
moment  and  HOMO-LUMO  energy  gap  of  the  optimized  isomeric  configurations  are 
listed  in  Table  10. 


31 


A:  Approved  for  public  release  (PA);  distribution  unlimited.  PA  Case  No:  88ABW-2012-0206  Date  Cleared:  19  November 
2014 


(a)  i-cooH 


%4 

r 

^  T' 


V  .23 

^*1^1 

/  ■' 

-^11 


*•, 

-4& 


'tf  « 

iU  ^ 


%' 


(b) 


Il-Indole 


tv 


_■# 


v 


.^Ig4 

'V  ■# 


Figure  23.  The  ground  state  configurations  of  Trp-Au32  complex. 
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Table  10.  The  calculated  structural  and  electronic  properties  of  the  isomeric 
configurations  of  the  Trp-Au32  complex. 


Isomers 

l-COOH  group 

W-Indole 

III-NH2 

Binding  energy/eV 

-2.72 

-1.86 

-1.8 

R(Cluster-Trp) 

- 

- 

- 

Au-N 

3.0 

3.31 

2.60 

Au-0 

2.58 

- 

2.9 

Au  —  C  /ndo/e-benzene 

2.48 

2.35 

- 

Au  —  C  /ndo/e-pyrrole 

2.38 

2.95 

- 

Au  -  H 

2.57 

2.71 

2.40 

Dipole  moment/D 

5.97 

3.78 

7.44 

HOMO-LUMO 
gap/e  V 

1.2 

1.51 

1.45 

Mulliken  charge/e 

Au  =  -0.062 

Trp  =  +0.062 

Au  =  -0.08 

Trp  =  +0.08 

Au  =  -0.13 

Trp  =  +0.13 

The  binding  energy  of  the  Trp-Au32Ciuster  compiex  is  calcuiated  as  the  difference  of 
the  energy  associated  with  the  equiiibrium  configuration  and  the  energies  associated 
with  the  isolated  ciuster  and  Trp.  As  noted  from  Table  10,  the  salt  bridge  structure 
invoiving  the  -COOH  group  appears  to  be  energeticaliy  more  favorable  than  the  two 
comparably  interacting  -indole  and  -NH2  groups.  The  larger  binding  energy  of  Trp-Au32 
in  -COOH  group  conjugation  appears  to  resuit  from  the  interaction  of  Au  with  carbonyl 
oxygen  atoms  (Figure  23a).  Besides  the  Ojrp-Au  interaction,  there  exists  a  weak  non- 
conventionai  Oirp-H-’-Au  interaction,  which  seems  to  provide  additional  stability  to  the 
Trp-Au  complex.  Figure  23a  aiso  shows  Au  atoms  interacting  with  /nc/ony/ functional 
groups,  including  both  of  its  pyrrole  and  benzene  rings  in  Trp.  There  are  some  structural 
changes  observed  for  Trp.  The  angle  C1-C2-C4  is  found  to  be  117-  compare  to  the 
initial  value  of  113.77-  and  the  angle  C2-C4-C9  becomes  111.9-  compare  to  its  initial 
value  of  105.45-.  However,  none  of  the  Hjrp  is  found  to  be  perturbed.  With  regards  to 
the  -indole  group  interaction,  Au  atoms  are  found  to  be  weii  interacting  with  indole 
group  of  Trp  (Figure  2b).  However,  no  noticeable  structural  changes  occurred  from 
indole  group  interaction.  The  angles  C1-C2-C4  and  C2-C4-C9  show  minor  aiterations 
and  are  found  to  be  at  114.54-  and  106.93-  respectively.  Unlike  the  -COOH  site 
interaction  having  an  -indole  contribution,  the  -NH2  site  interacts  with  the  nitrogen  of 
the  amine,  oxygen  of  the  carboxyl  group  and  hydrogen  of  the  carbon  backbone  (Figure 
23c).  Figure  2c  shows  the  double  coordination  of  Au  with  -COOH  and  -NH2  groups  and 
the  angles  C1-C2-C4  and  C2-C4-C9  appear  at  118.76-  and  112.69-  respectively. 
Additionally,  the  angle  between  the  -COOH  group  and  -NH2  group  in  Trp  was  reduced 
to  106.33-  from  its  initiai  value  of  109. 52-. 

In  protein  chemistry,  structures  invoiving  the  -COO-  group  of  an  amino  acid  are  known 
as  sait  bridge  structures,  whereas  interactions  invoiving  the  -indole  ring  are  represented 
in  charge  solvent  structures  [33].  Following  this  nomenclature,  the  Trp-Au32  isomer 
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formed  through  the  -COOH  group  interaction  can  be  ciassified  as  a  sait  bridge  structure 
but  as  noted  from  Figure  23a  and  Table  10,  the  Au  cluster  is  equally  interacting  with  the 
-indole  group-.  Hence,  the  isomeric  configuration  involving  the  -COOH  group  may  be 
considered  to  have  characteristics  of  both  a  salt  bridge  as  weii  as  a  charge  soivation 
structure,  showing  maximum  binding  energy  (Table  10).  Likewise,  the  isomeric 
configuration  invoiving  the  -indole  group  may  be  considered  as  a  charge  soivation 
structure.  The  -NH2  group  isomer  shows  double  coordination  with  the  -COOH  and 
-NH2  groups  (Figure  23c).  Thus,  the  caiculated  results  predict  the  ground  state 
complex,  which  is  in  agreement  with  the  experimentai  observation  with  a  preference  for 
a  mixed  sait  bridge  and  a  charge  soivation  type  bonding  interaction  [23,24]. 

The  Muiiiken  charge  anaiysis  of  the  isomeric  configurations  of  the  Trp-Au32  compiex 
reveais  a  small  charge  transfer  during  conjugation  and  is  noted  in  theTrp  to  Au  cluster 
at  all  three  sites  of  interaction.  The  results  therefore  excluded  a  charge  transfer  as  the 
main  mechanism  for  the  Trp-Au32  conjugation  for  the  -COOH  and  -indole  group 
interaction.  However,  the  charge  transfer  is  comparativeiy  higher  in  the  case  of  the  -NH2 
group  interaction  that  may  govern  the  conjugation.  This  is  in  agreement  with  the 
reiatively  iarger  dipole  moment  of  7.44  debye  (D)  for  this  complex  during  the  -NH2 
group  interaction  compared  to  the  native  Trp,  suggesting  the  rearrangement  of  the 
electron  density  when  the  complex  forms.  The  comparatively  large  dipole  moment  of  the 
-COOH  group  compared  to  that  of  the  -indole  group  is  co-related  with  a  iarger 
structurai  rearrangement  evidenced  during  the  formation  of  -COOH  group  conjugation; 
whereas,  no  such  structurai  perturbation  is  noted  for  the  -indole  group  conjugation. 
Interestingly,  both  the  -COOH  and  -indole  interaction  sites,  showing  nearly  the  same 
values  of  small  charge  transfers,  share  a  common  factor  of  an  indole  group  during 
conjugation  with  AU32  (Figure  23).  The  reiatively  stronger-than-expected  interactions 
from  the  indole  group  in  the  computed  results  occur  as  a  result  of  the  indole  group 
donating  a  tt  electron  cloud  to  the  Au  atoms.  The  planarity  of  indole  also  maximizes  the 
deiocalization  interaction.  These  effects  were  observed  in  previous  studies  [34].  Thus, 
our  calculation  suggests  that  the  preferred  nature  of  the  bonding  was  dominated  by  tt 
interactions  with  only  a  small  charge  transfer. 

Figure  24  shows  the  charge  density  piots  projected  aiong  the  001  piane  for  aii  the 
isomeric  configurations  of  the  Trp-Au32  compiex  (Figure  24b)  compare  to  native  Trp 
and  AU32  (Figure  24a).  Figure  24b  ciearly  shows  that  the  interaction  and  charge 
redistribution  are  taking  place  with  the  indole  ring  and  carbonyi  oxygen  in  Trp-Au32 
complex  for  I  -COOH  group  interaction,  corroborating  previous  resuits.  For  the  W- indole 
group  compiex  in  Figure  24b,  the  charge  density  plot  shows  a  weak  charge 
delocalization  with  the  -indole  group,  indicating  existing  tt  electronic  interactions  for 
Trp-Au32  conjugation.  The  charge  density  piot  for  III-NH2  in  Figure  24b  reveals  an 
electronic  charge  redistribution  during  -NH2  group  interaction  that  is  higher  than  the 
other  two  mentioned  interactions.  Hence,  charge  density  results  provide  evidence  for 
the  nature  of  interactions  as  predicted  in  Figure  24. 
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(a)  Total  charge  density  of  native  AU32  and  Trp 


AU32 


Trp 


(b)  Total  charge  density  of  Trp-Au32  complex 


Figure  24.  Total  charge  density  plot  of  AU32,  Trp  and  isomeric  configurations  of  Trp-Au32 
complex  projected  along  001  plane.  A  superimposed  ball  and  stick  model  identifies  the 
atoms  in  the  complex.  Total  charge  density  is  shown  at  1/3rd  of  its  maximum  isosurface 
value. 


The  previous  experimental  study  identifies  the  binding  site  of  Au  nanoparticles  in  BSA 
containing  Trp  in  near  vicinity,  resulting  in  fluorescence  quenching  [23,24].  Since  the 
spectral  signatures  of  the  protein-Au  conjugate  are  dominated  by  the  tt-tt*  transition  of 
Trp’s  indole  functional  group,  the  formation  of  the  complex  could  modify  the  electronic 
states  or  introduce  new  states  in  Trp,  resulting  in  loss  of  fluorescence  energy.  In  order 
to  examine  this,  we  plotted  the  highest  occupied  molecular  orbital  (HOMO)  and  lowest 
unoccupied  molecular  orbital  (LUMO)of  the  pristine  Trp,  AU32  and  Trp-Au32  isomers  in 
Figure  25. 
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Figure  25.  The  contour  plots  of  molecular  orbitals  of  pristine  Trp,  AU32  and  Trp-Au32 
conjugate  isomers  shown  at  1/2nd  of  its  maximum  isosurface  value. 

As  seen  in  the  figure,  the  HOMO  and  LUMO  of  pristine  Trp  are  localized  at  the  indole 
functional  group.  However,  significant  modification  in  the  nature  of  molecular  orbitals, 
and  hence  the  charge  re-distribution  happens  upon  the  Trp-Au32  conjugate  formation. 
The  orbitals  associated  with  Au  atoms  appear  as  low-lying  LUMO  in  the  energy 
spectrum  of  all  the  three  isomers.  This  LUMO  may  potentially  act  as  an  acceptor  from 
the  excited  fluorescence  state  of  Trp  in  the  complex,  thus  reducing  the  fluorescence 
intensity.  The  higher  LUMOs  i.e.  LUMO-h1,  LUMO-h2  are  delocalized  over  Au  atoms  as 
well  as  on  the  interacting  Trp  atoms.  The  delocalized  higher  LUMOs  indicate  that  the 
hybridization  in  Trp-Au32  yields  orbitals  of  a  mixed  metal-molecule  character. 

Appearances  of  Au  localized  LUMO,  hybrid  LUMO-i-1  and  LUMO-i-2  orbitals  are  further 
confirmed  by  the  total  density  of  states  (DOS)  of  isomers  as  shown  in  Figure  26.  The 
DOS  of  Trp-Au32  isomers  show  hybrid  states  near  to  Fermi  energy  that  are  introduced 
in  the  system  by  Au,  which  most  likely  causes  fluorescence  quenching  as  measured  in 
our  experiments. 
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Figure  26.  Density  of  states  of  Trp,  AU32  and  Trp-Au32Compiex  isomers:  (a)  l-COOH, 
(b)  ll-indoie  and  (c)  III-NH2.  The  HOMO  is  aligned  to  origin  in  energy  scale.  The  low 
lying  hybrid  excited  states  introduced  by  Au  are  marked  in  pink  circle. 


The  aromatic  organic  molecules  are  reported  to  show  signatures  of  weak  chemisorption 
at  gold  surfaces,  producing  the  bonding-  and  anti-bonding-llke  hybrid  orbitals  below  and 
above  the  middle  of  the  metal  d  band,  as  Illustrated  in  the  Newns-Anderson  model 
(Figure  27f)  [35].  The  model  shows  that  the  metal  d  electron  band  Interacts  with  highest 
occupied  electron  level  of  the  interacting  atom  and  produces  bonding  and  anti-bonding 
levels.  Since  the  chemisorption  is  weak,  the  signatures  for  these  levels  are  also 
expected  to  be  weak.  Thus,  the  nature  of  hybridization  at  the  interacting  interface  of 
Trp-Au32  isomers  is  determined  in  detail  from  the  computed  projected  density  of  states 
(PDOS). 

Figure  27  represents  projected  density  of  states  (PDOS)  centered  on  atoms  involved  in 
the  hybridization  during  Trp-Au32  complex  formation,  compared  to  their  native  state. 
Figure  27a  shows  PDOS  of  native  Au  d  electrons  before  conjugation,  showing  strong 
peaks  in  the  -6  to  -2  eV  energy  range  and  is  apparently  localized.  Figure  27b  shows  the 
p  orbitals  on  native  Trp  atoms  participating  in  the  hybridization.  The  PDOS  of  carbon 
atoms  in  the  indole  ring  are  summed  up  and  presented  together  as  indole  C,  while  the 
indole  nitrogen  N1  is  shown  separately.  The  structure  of  Trp  with  marked  atom  numbers 
is  shown  in  Figure  22b.  The  indole  carbon  atoms  show  sharp  peaks  due  to  the 
presence  of  a  delocalized  tt  electron  cloud  in  the  aromatic  ring.  The  nitrogen  atom  N1  is 
also  a  part  of  this  aromatic  tt  electron  system  and  therefore  found  to  follow  the  same 
trend  of  the  indole  C  spectrum  with  less  intensity  (Figure  27b).  There  is  a  strong  peak 
noted  around  -1 .2  eV,  which  is  attributed  to  the  lone  pair  of  the  nitrogen  atom  N2  of  the 
amine  group  in  Trp.  The  oxygen  01  also  appears  with  a  strong  peak  around  -3.5  eV, 
related  to  Its  lone  pairs.  The  oxygen  atom  02  appears  with  two  strong  peaks  around  - 
3.5  and  -2  eV  related  to  its  lone  pair  and  carbonyl  tt  electrons  respectively. 

Figure  27c  represents  the  PDOS  of  hybridizing  atoms  of  Trp-Au32  complex  at  -COOFI 
site.  Following  the  Newns-Anderson  model,  the  energy  range  of  -8  to  -4  eV  is  assigned 
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to  the  bonding  region  and  -2  to  0  eV  is  assigned  to  the  anti-bonding  region.  The  peaks 
corresponding  to  Au,  oxygen  01,  oxygen  02,  indole  C  and  N1  have  been  modified 
compared  to  the  initial  one  in  Figure  27b,  confirming  their  participation  in  hybridization 
and  their  contribution  to  bonding  and  anti-bonding  orbitals.  Similar  to  this.  Figure  27d 
represents  the  PDOS  of  the  Trp-Au32  complex  for  the  -indole  interacting  site,  showing 
participation  of  Au,  C2,  indole  C  and  N1  in  the  hybridization  and  formation  of  bonding 
and  anti-bonding  peaks  from  -8  to  -4  eV  and  -2  to  0  eV  respectively.  Figure  27e  shows 
the  PDOS  of  hybridizing  atoms  of  the  Trp-Au32  complex  during  the  -NFI2  site 
interaction.  The  lone  pair  of  N2  has  been  largely  modified  as  would  be  expected  for  this 
interaction  as  with  02.  The  01  also  seems  to  participate  in  this  interaction,  as  its  orbital 
signatures  have  been  modified.  The  anti-bonding  region  of  the  hybridizing  system  in  all 
three  interacting  site  is  actually  crossed  by  the  fermi  level;  therefore  the  region  has  been 
found  to  be  less  populated.  The  presence  of  the  underlying  broad  metal  s  band  affects 
the  nature  of  hybrid  bonding  and  anti-bonding  orbitals  and  generates  spread  peaks 
rather  than  single  localized  peaks. 

It  is  clearly  demonstrated  in  the  PDOS  plots  of  Trp-Au  isomers,  that  the  peak 
corresponding  to  the  LUMO  (~1eV)  of  the  isomer  is  mainly  centered  at  the  interacting 
Au  atoms,  though  the  higher  excited  levels  show  a  mixed  Au-Trp  character  as  evident 
by  molecular  orbital  (Figure  25)  and  DOS  analyses  (Figure  27). 
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Figure  27.  The  projected  density  of  states  (PDOS)  of  interacting  atoms  in  (a)  native 
Au32,  (b)  native  Trp,  (c)  -COOH  group,  (d)  -indole  group,  (e)  -NH2  group  isomers  of 
Trp-Au32  and  the  (f)  Newns- Anderson  bonding  scheme.  Fermi  levels  are  aligned  at  the 
origin  of  the  energy  scale 
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Conclusion 


First  principies  caiculations  based  on  DFT  were  performed  on  Trp,  AU32  and  Trp-Au32  in 
order  to  study  the  details  of  the  nature  of  the  interaction  between  tryptophan  and  goid. 
We  find  that  preference  of  the  binding  involving  both  carboxyi  and  indole  functionai 
groups  yielding  the  mixed  sait  bridge  and  charge  soivation  structure  to  be  energetically 
more  favorable  and  mediated  by  stronger-than-expected  tt  interactions  offered  by  the 
indole  group.  The  present  study  does  support  the  experimentai  finding  of  a  ground-state 
Trp-Au  complex  formation.  However,  the  present  study  depicts  tt  interactions  along 
with  smail  ionic  interactions,  instead  of  van-der  Waals  between  Trp  and  Au 
nanoparticies.  The  hybridization  during  the  formation  of  the  stable  Trp-Au  complex  is 
anticipated  to  modify  the  nature  of  Trp’s  moiecular  orbitals  by  introducing  hybrid  states 
associated  with  Au  as  iow  lying  excited  states.  The  PDOS  anaiysis  suggests  that  these 
hybrid  orbitais,  produced  by  the  coupiing  of  the  narrow  Au  d  band  with  Trp  atomic 
orbitais  have  a  mixed  metal-moiecuie  character  and  are  delocalized  at  the  interface. 
The  non-radiative  energy  transfer  to  these  states  is  likely  to  cause  the  experimentally 
observed  fluorescence  quenching  in  Trp. 
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2.  Lipid-Nanoparticle  Interactions:  DMTAP,  DMOC  and  DMPE  Bilayers 
Introduction 

A  biological  membrane  separates  a  cell  from  its  surrounding  environment,  and  its 
function  depends  on  the  continuity  and  integrity  of  the  bilayer  structure,  since  over  70% 
weight  of  membrane  belong  to  a  lipid  bilayer  [36].  Consequently,  knowledge  of  the 
structure  of  a  lipid  bilayer  is  the  key  element  to  understand  the  functionality  of  a 
biological  membrane,  especially  the  tertiary  structure  which  can  be  obtained  by  the  X- 
ray  diffraction  measurements  [37-43]. 

In  recent  years  theoretical  methods  have  become  one  of  the  important  tools  to 
investigate  the  physical  properties  of  a  biological  system.  These  methods  are  generally 
based  on  molecular  dynamics  (MD)  and  Monte  Carlo  (MC)  methods,  and  can  provide 
three-dimensional  real-time  imaging  of  the  complex  biological  system  and  simulate 
many  important  phenomena  including  self-assembly,  phase  transition  [43-49],  and  lipid- 
molecule  interactions  [50-55],  and  the  lipid  bilayer  system  [56-62].  We  note  that  these 
atomistic  simulation  methods  are  quite  limited  when  it  comes  to  obtaining  the  electronic 
structure  of  lipid  bilayers.  Furthermore,  their  accuracy  depends  on  the  description  of  the 
force  field  parameters  obtained  by  the  empirical  methods. 

In  this  study,  we  employ  an  alternative  way  of  calculating  both  structural  and  electronic 
properties  of  lipid  bilayers  which  does  not  require  the  knowledge  of  the  force  field 
parameters.  Our  first-principles  method  is  based  on  density  functional  theory.  Recently, 
Mashaghi  and  co-workers  used  first  principles  method  to  investigate  hydration  affection 
on  molecular  and  electronic  structure  of  phospholipid  [63].  It  is  to  be  noted  here  that  the 
overall  aim  of  this  project  is  to  investigate  interactions  between  nanoparticles  and  cell 
membranes,  which  are  central  to  understanding  the  potential  deleterious  effects  of 
nanoparticle  on  living  cells. 


Computational  method 

Since  calculations  based  on  first  principles  methods  are  computationally  intensive,  our 
challenge  will  be  how  to  simplify  the  lipid  bilayer  structure  for  simulation  without  losing 
its  essential  structural  features  which  are  important  in  the  understanding  of  its 
functioning.  Assuming  a  2-dimensional  periodic  structure  of  the  lipid  layer,  we  employed 
the  framework  of  the  generalized  gradient  approximation  (GGA)  to  density  functional 
theory  using  the  Vienna  ab  initio  Simulation  Package  (VASP).  The  projector 
augmented-wave  (PAW)  [64-67]  potentials  and  plane  wave  basis  sets  together  with  the 
exchange  and  correlation  functional  forms  proposed  by  Perdew  and  Zunger  [69,70] 
were  used.  During  the  optimization  procedure,  the  size  of  the  super  cell  as  well  as  the 
internal  coordinates  of  all  the  atoms  was  relaxed.  The  energy  cutoff  of  the  plane  wave 
was  taken  to  be  400  eV. 


Results  and  Discussion 
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DMTAP  (Figures  28  and  29)  belongs  to  a  class  of  double-chained  cationic  lipid 
amphiphiles,  and  has  shown  considerable  promise  as  potentially  effective  DNA 
transfection  agents.  The  lipid  structure  has,  in  general,  two  types  of  chemical  groups. 
One  chemical  group  brings  the  positive  charge  such  as  the  choline,  amine  etc.,  whereas 
the  negative  charge  is  associated  with  the  phosphate  group.  All  of  the  positive  charges 
in  the  lipid  layer  lie  on  a  plane  or  shell,  and  all  of  the  negative  charges  lie  in  the  other 
plane  or  shell.  These  two  charge  plane  or  shells  make  up  the  electric  double  layer  which 
mainly  determines  the  functionality  of  a  cell  membrane. 


Figure  28.  A  ball  and  stick  model  for  DMTAP  (O  :  red,  C :  gray ,  N:  yellow,  H :  blue) 


Figure  29.  A  ball  and  stick  model  for  DMTAP  layer.  (O :  red,  C :  gray ,  N:  yellow,  H : 
blue) 


In  our  simulation  model,  the  super  cell  of  the  DMTAP  lipid  layer  has  a  triclinic  Bravais 
Lattice.  The  average  area  per  lipid  molecule  is  calculated  as  follows: 


A  =  b  X  c  X  sin  a 


The  calculated  values  of  the  parameters  for  the  optimized  DMTAP  unit  cell  and  the 
experimental  results  [70]  are  shown  in  Table  1 1 . 
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Table  1 1.  The  calculated  structural  properties  of  DMTAP. 


b(A) 

c(A) 

a(°) 

Area/lipid  (A2) 

Volume  (A3) 

This  work 

4.4 

8.0 

96.2 

34.9 

904 

Experiment  [36] 

5.1 

7.3 

92 

37 

917 

We  note  that  the  average  area  per  lipid  is  one  of  most  fundamental  characteristics  of 
lipid  bilayers.  The  excellent  agreement  between  our  calculated  results  with  the 
experimental  value  successfully  validates  the  modeling  elements  employed. 

The  lipid  bilayers,  phosphatidylcholine  (PC)  and  phosphatidylethanolamine  (PE), 
represent  the  abundant  phospholipids  class  in  animal  and  plant  tissues,  and  can  be  the 
major  lipid  class  in  microorganisms.  In  the  present  study,  we  take  1 ,2-Dlmyrlstoyl-sn- 
glycero-3-phosphorylcholine  (DMPC)  and  1 ,2-Dimyristoyl-sn-glycero-3- 
phosphorylethanol-  amine  (DMPE)  as  the  representative  members  of  the  PC  and  PE 
types  shown  in  Figures  30  and  31 .  Note  that  the  difference  between  the  PE  and  PC  is 
the  N  atom  in  the  hydrophilic  head  -  it  is  connected  to  three  methyl  groups  in  DMPC, 
but  is  connected  to  three  H  atoms  in  DMPE. 


Figure  30.  A  ball  and  stick  model  for  DMPC  layer.  (O  :  red,  C  :  gray ,  N:  yellow,  H :  blue) 
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Figure  31.  A  ball  and  stick  model  for  DM  PE  layer.  (O :  red,  C :  gray ,  N:  yellow,  H :  blue) 


In  order  to  make  calculations  computationally  feasible,  we  first  performed  calculations  to 
investigate  the  number  of  C  atoms  required  to  represent  the  “tail”  of  a  fatty  acid,  which 
was  a  long  hydrocarbon  chain.  Table  12  lists  the  results  of  such  calculations  where  we 
find  that  the  area  associated  with  8  C  atoms  representing  the  "tail"  of  the  fatty  acid  was 
reasonably  enough  for  further  simulation  of  DMPC  and  DMPE.  This  was  consistent  with 
the  experimental  study  which  reported  the  average  area  per  lipid  molecule  remains 
nearly  the  same  for  the  long  hydrocarbon  chain. 


Table  12.  The  average  area  per  lipid  molecule  as  a  function  of  number  of  C  atoms 
representing  the  "tail"  of  a  fatty  acid.  Units  are  in  A^. 


C  atoms 

representing  the  "tail" 

DMPC  ^ 
Area/lipid  (A^) 

DMPE 

Area  (A^) 

2 

33.60 

32.49 

4 

36.97 

33.88 

8 

37.43 

35.64 

10 

37.62 

36.50 

14 

37.60 

39.70 

Table  13.  The  calculated  structural  properties  of  DMPC  and  DMPE. 


b(A) 

c(A) 

a(°) 

Area/lipid(A^) 

Volume  (A^) 

DMPC:  This  work 

4.4 

8.6 

95.1 

37.6 

954.1 

DMPC:  Experiment[37] 

4.71 

8.46 

90.0 

39.8 

- 

DMPE:  This  work 

4.29 

9.26 

93.1 

36.5 

1154.9 

DMPE  Experiment  [38] 

3.73 

9.70 

90.0 

36.1 

- 

Both  the  DMPC  and  DMPE  can  be  represented  by  the  triclinic  Bravais  lattice  whose 
parameters  are  given  in  Table  13.  Considering  that  the  area  per  lipid  becomes  smaller 
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as  the  temperature  decreases  [72,  73],  the  calculated  value  at  OK  follows  the  trend  as 
the  experimental  value  was  obtained  at  283K.  A  comparison  between  the  calculated 
and  experimental  values  of  the  lattice  parameters  shows  a  good  agreement  for  the  lipid 
layers  considered.  We  note  that  the  previous  calculations  using  the  MD  method 
generally  found  reported  the  per-area  to  be  50-60  for  DMPC  and  DMPE  lipid  bilayers 
at  the  low  temperature  [71]. 

It  is  well  known  that  many  of  the  physical  and  chemical  properties  of  a  lipid  bilayer  can 
be  determined  via  knowledge  of  the  electron  density.  Figures  32  and  33  show  the 
electron  density  contour  plots  of  the  electric  double  layer  part  of  DMPC  and  DMPE, 
respectively. 


Figure  32.  The  eiectron  density  contour  plots  of  DMPC. 
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Figure  33.  The  electron  density  contour  plots  of  DMPE. 


In  these  systems,  the  contour  plots  show  that  P  has  lost  the  electron  density  to  nearby 
oxygens.  On  the  contrary,  N  is  more  electronegative  than  C,  thus  the  region  near  N 
atoms  shows  excess  in  the  electron  density. 

In  order  to  mimic  the  cell  live  environment,  we  now  add  the  water  molecules  in  our 
simulation  model  for  DMPC  and  DMPE.  The  corresponding  optimized  structures  of  lipid 
monolayers  are  shown  in  Figures  34  and  35. 


Figure  34.  A  ball  and  stick  model  for  the  DMPC:  H2O  system. 


(O:  red,  C:  green,  P:  yellow,  N:  blue  H:  golden) 
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Figure  35.  A  ball  and  stick  model  for  the  DMPC:  H2O  system. 
(O:  red,  C:  green,  P:  yellow,  N:  blue  H:  golden) 


For  a  lipid  biiayer,  the  magnitude  of  the  eiectric  dipole  moment  mainly  determines  the 
degree  of  activity  of  the  membrane  [74,  75].  For  DMPC  and  DMPE  systems,  the 
calculated  values  are  tabulated  in  Tabie  14  considering  the  center  of  ceii  to  be  (0.5,  0.5, 
0.5)  as  the  reference  of  point. 


Table  14.  The  calculated  dipole  moments  of  DMPC  and  DMPE  monolayers. 


a  (exA) 

(3  (exA) 

Y  (ex A) 

Energy  (eV) 

DMPC 

isolated 

hydrated 

-0.851 

-0.842 

0.115 

0.123 

-2.561 

-2.544 

0.0286 

0.0283 

DMPE 

isolated 

hydrated 

-0.082 

-0.411 

0.011 

0.015 

-1.032 

-3.184 

0.0042 

0.0424 

For  the  bilayer  system,  our  calculated  values  will  be  -1 .32  and  -0.78  kcal/mol  for  DMPC 
and  DMPE,  respectively.  A  theoretical  study  using  the  iinearized  Poisson-Boltzmann 
equation  reported  the  vaiue  of  the  eiectrostatic  potential  at  the  bilayer  mid-plane  to  be  - 
1.21,  for  DMPC  bilayer  [74,  75].  Furthermore,  the  results  find  the  presence  of  the  water 
moiecules  to  be  the  determining  factor  in  stabiiizing  the  functioning  structure  of  the 
DMPE  bilayer,  and  facilitates  a  significantly  large  change  in  the  dipoie  moment  along 
the  Y  direction.  The  resuits  therefore  predict  a  reiatively  high  activity  of  the  hydrated  iipid 
layer  of  DMPE  relative  to  that  of  DMPC  for  nanoparticles. 


Conclusion 

Considering  that  understanding  of  the  interactions  between  nanoparticles  and  a 
biomembrane  is  centrai  to  predict  the  potential  deleterious  effects  of  nanoparticie  on 
living  ceiis,  we  performed  simuiation  of  iipid  biiayers  using  first  principles  methods.  The 

48 

A:  Approved  for  public  release  (PA);  distribution  unlimited.  PA  Case  No:  88ABW-2012-0206  Date  Cleared:  19 
November  2014 


calculated  results  show,  for  the  first  time,  the  accuracy  and  reiiabiiity  of  such  methods  in 
calcuiating  the  structural  property  of  for  DMTAP,  DMPC  and  DMPE  bilayer  systems.  We 
also  examine  the  solvation  effects,  and  find  that  the  presence  of  the  water  moiecuies  is 
essentiai  in  stabilizing  the  configuration  of  the  DMPE  bilayer.  It  also  facilitates  a 
significantiy  large  change  in  its  dipole  moment  along  the  y  direction.  The  results 
therefore  predict  a  reiativeiy  high  activity  of  the  hydrated  iipid  iayer  of  DMPE  reiative  to 
that  of  DMPC  for  nanoparticies. 
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